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. We present exact analytical zero-energy solutions for a class of smooth decaying potentials, show- 

^N) ' ing that the full confinement of charge carriers in electrostatic potentials in graphene quantum dots 

' and rings is indeed possible without recourse to magnetic fields. These exact solutions allow us 

O , to draw conclusions on the general requirements for the potential to support fully confined states, 

■ including a critical value of the potential strength and spatial extent. 

0> ' PACS numbers: 73.22.Pr, 73. 21. La, 03.65.Ge, 03.65.Pm 

(N : 

I. INTRODUCTION 

^ : 

. There is a widespread belief that electrostatic confinement of charge carriers in graphene is not possible due to the 
c/5 ' efi^ect of Klein tunneling,! by which carriers inside a potential well may couple to states outside the potential via empty 
^ ', states in the hole band. The conventional interpretation of this effect in graphene^''^ is that the conservation of chirality 
^ ■ forbids backscattering for normally incident particles, therefore the transmission probability is unity irrespective of 
[ the barrier height. However, the effect is diminished for particles not normally incident on the barrier which can be 
^ • achieved by particles possessing a finite transversal wave vector^ 

' As a result of Klein tunneling, many considerations of axially symmetric potential wells maintain that confinement 
inside electrostatic quantum dots is impossible;^ Instead, focus has been placed on formation of confined states within 
graphene by the application of magnetic fields perpendicular to the graphene plane^""— or on engineering specialized 
devices which introduce mass-like terms.-WrJ^ Nonetheless, previous authors have studied the manipulation of charge 
^ ' carriers by electrostatic fields such as the creation of quasibound states in quantum dots, whose lifetimes are long 
but not infinite, as certain angular momenta correspond to trajectories incident on the potential barrier with low 

■ transmission probabilities i^ii^ An alternative interpretation of long-lived states in electrostatic dots is that a wcU- 
^ ' tuned structure can generate wavefunctions outside the dot which interfere destructively 

It is well known that an ideal two-dimension system with linear energy dispersion possesses a density of states 
' that vanishes in the limit of zero energyji^ This raises the possibility of confining states by careful manipulation of 
OO , the Fermi energy by application of a back-gate voltage, allowing the energy of the system to coincide with the Dirac 
points. In this scenario, no states exist outside of an electrostatic potential well for bound states to tunnel into. This 
I approach has previously been developed with reference to electron waveguidesj^ for tunneling selection of charge 

■ carriers through an n-p junction,— and for sharply terminated wells of constant depth.— In this paper we look at 
full confinement via this method in smooth potentials with circular symmetry. We would like to emphasize that our 

" ^ I I analytic results are obtained within an idealized model neglecting impurity- and substrate-induced disorder, local 
^ ' strain distribution, ripples, and trigonal warping as well as many-body and finite-size effects, which all influence the 
energy spectrum and transport properties of realistic graphene flakes, especially near the charge neutrality pointji^ 
The application of an external fleld with circular symmetry would be possible by close proximity of a scanning 
I tunneling microscope (STM) tip. The presence of bound modes can then be probed by performing conductivity 
measurements, in which electrons may propagate between terminals via the states created by the potential. Similar 
experiments relevant to the creation of localized states have recently been performedJ^i^ 

Quantum dots are important for the understanding of fundamental physics and also for applications, for example, 
being crucial to progress in spintronics and quantum computation. Applications of graphene to these disciplines is 
doubly important: flrst because graphene is considered to be a system where zero-energy Majorana fermions may 
appear,-! g^j^^j second because quantum dots in graphene have potential for use as spin qubits.^^ Notably, most of 
the results of this work are also applicable to the Dirac-like states on the surface of topological insulators, ^■^ which is 
another exciting area of contemporary condensed matter physics. 

The Hamiltonian operator in the massless Dirac- Weyl model for graphene, which describes the motion of a single 
electron in the presence of an external axially symmetric potential U{r), is2ii2^ 
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H = vpcr ■ p + U{r), 
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where vp « c/300 is the Fermi velocity of the charge carriers, <t = {<7x,<Jy) are the Pauh spin matrices, and p = —ihSJ 
are the hnear momentum operators. We consider smooth confining potentials, which do not mix states in the two 
nonequivalent graphene valleys, and deal with one valley only. All our results can be easily reproduced for the other 
valley. For axial symmetry, we transform into cylindrical coordinates (r, 0) and separate the variables by using the 
following ansatz for the two-component wavefunction 



2ti V XB{r) 



where m is the integer- valued angular momentum number and the subscripts A and B refer to the sublattices of the 
graphene honeycomb lattice. This choice of wavefunction leads to two coupled first order differential equations for 
the radial wavefunction components xa.b 



d m + 1 
dr r 
d m 
dr r 



iXB = {e-V{r))xA, 

(3) 

iXA = {£-V{r))xB, 



with V{r) = U{r)/hvF and e = E/hvp, where E is the energy eigenvalue. It is most convenient to solve these 
equations by formulating a second order equation for one component and reusing one of Eqs. ^ to find the other 
wavefunction component. It is worth noting that the radial wavefunction components transform into each other on 
replacing m — —(to -|- 1), which is important when considering special values of to. 

The rest of the paper is organized as follows. In SecUll we look into the long-range behavior of zero-energy states in 
power-law potentials decaying at large distances faster than the Coulomb potential and show that only wavefunctions 
corresponding to a non-zero angular momentum (vortices) are square-integrable across the whole graphene plane. 
Section [IIII contains the main results of our work - hitherto unknown analytic solutions for fully confined zero-energy 
states in a class of circularly symmetric decaying potentials. In Sec. II VI we discuss the relevance of our results to STM 
of graphene. Finally, in Sec.|V]we summarize our results and formulate possible further applications and developments 
of our work. 



II. LONG-RANGE BEHAVIOR IN A DECAYING POTENTIAL 



The Coulomb problem for Dirac fermions (the Dirac-Kepler problem) in two dimensions has been studied for some 



time, - and there has been a resurgence of interest in this problem since the discovery of graphene. 



Despite its 



apparent beauty, the ideal Coulomb problem has somewhat questionable relevance to the reality of graphene since 
the impurity potential should have a short-range cutoff of the Ohno typc"^" as well as a faster than 1/r long-range 
decay due to either screening or an image charge appearing in the metallic back gate, which is necessarily present in 
graphene-based devices. 

For a more realistic quickly decaying potential when a system is fixed by the back-gate at the Dirac point energy 
(_E = 0), there is no coupling to the continuum and, as we show, a single fully bound state might appear. General 
properties of the fully confined solutions can be understood from their long-range behavior. Thus we consider the 
general potential given by 



Vir) = Vor- 



(4) 



where the rate at which the potential falls off is characterized by p > 1. For this potential the differential equation 
obeyed by the upper wavefunction component is 



Xa 



l+p 
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XA 



0. 



(5) 



This equation can be recast as the Bessel differential equation. In general the solution is a linear combination of 
Bessel functions in the form 



XA{r) = r 



-p/2 



Vr 



ClJa ( — P ) + C2J-a , 

p-1 J 



,.1-? 
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where a = (p + 2to) /(2 — 2p) is the order of the function. Then it follows from Eqs. ([3]) that the lower radial 
wavefunction component is 



XB{r) = ir 



= ,>-p/2 



p-l 



(7) 
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Now we look into the long-distance behavior of these solutions in order to gain insight into the behavior of carriers 
in a general confining potential. For large radial coordinate, the variable of the Bessel functions tends to 0, thus the 
desired behavior is dominated by the first-order term in the Maclaurin expansion. For m > the asymptotic solution 
becomes 

- ^B{r) ^ (8) 

Similarly for m < — 1 the asymptotic solution is 

X^M^Cir™, XB(r)^*Ciri+"-f. (9) 

From this asymptotic behavior of the wavefunctions we can deduce that all states in any smoothly decaying potential 
are normalizable, with the exception of angular momentum states m = and m = —1, reflecting the m — > — (m -|- 1) 
symmetry. In these cases the integral of the probability distribution contains a logarithmically divergent term. All 
other angular momenta describe quasiclassical trajectories for which charge carriers are not normally incident on 
the confining potential (vortices), consistent with the idea that Klein tunneling in graphene is maximized at normal 
incidence and suppressed at other angles . ^'^'^'^^ 

It is worth noting that for p equal to an even integer, for at least a subset of m, the two solutions given in Eq. ^ are 
not linearly independent and the second independent solution for the wavefunction component is instead a Neumann 
function with the same order (a) as the first. However the asymptotic behavior of this function is the same as for the 
Bessel function with order —a, so the conclusions we draw are unchanged. 



III. SMOOTH QUANTUM DOTS AND RINGS 

We are interested in the presence of confined states when the system is fixed at the Dirac point energy, consistent 
with e = in the low-energy formulation. We consider the class of smoothly-varying radial potentials 

where Vb and d parametrize the strength and width of the potential respectively and k determines the sharpness 
of the ring. Most notable is the case where fc = 0, corresponding to a potential energy profile of the well-known 
Lorentzian shape. The benefit of a potential in this form is that it contains parameters that allow it to be fitted to 
experimentally realizable potentials and that it is both regular at the origin and short-range. The smoothness of the 
potential also allows one to neglect the intervalley scattering that occurs for tunneling problems, contrary to previous 
considerations of quantum dots with sharp boundariesj^'i^'i^ 

For simplicity, we first consider the case of the quantum dot, fc = 0. For the upper component of the wavefunction 
we seek a solution in the form 

XA{r) ^ Vir), (11) 

which leads to the second-order differential equation for 77 (r), 



r l + (r/d)2 J I"- >^ \ l + {r/dy^ 



V"ir) + ( + , ] V'ir) + { ) = 0, (12) 
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and changing to a new variable z — ir/d)^ we obtain an equation for 77(2), 



as 



z[l + zfr^"{z) + {l + z)[{m + l){l + z) + z\ri'{z)+{^^ 7,(2) = 0. (13) 
This differential equation has the following solution for tti > 0, which can be written using appropriate identities;^ 

77(z) = ci 2^1 in, -n; m + 1; ^ j , (14) 



where n = VQd/2 must take a non-negative integer value in order to terminate the hypergeometric function and Ci 
is a normalization factor. We can see that the solution is unchanged upon replacing Vq with — Vq, as required for 
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electron- hole symmetry. In order to have physically meaningful solutions we require that n > m, ensuring that the 
radial wavefunctions do not diverge. 
For m < — 1 the solution is 



r]{z) = C2 2F1 



1 



'1 + z 



(15) 



for which the corresponding requirement of the presence of physical wavefunctions is n > — m. 

The states formed by the presence of a Lorentzian-type electrostatic potential acting on a graphene sheet do indeed 
fulfill the expectation presented in Sec. [TTl that the states with m = and m — —1 are extended states with a 
nonintegrable probability density. Thus, the first normalizable state occurs for m = 1 (or m — —2), which requires 
n = 2 and so the threshold for fully confined states is Vod = 4. When the potential is fixed, charge carriers fill up 
modes enumerated by to. The apparent asymmetry between positive and negative values of to is removed when the 
other graphene valley is considered. 

A similar procedure allows one to solve the coupled equations with the potential describing electrostatic rings, 
namely, the potential ((T0| with k ^ 0. The general solution in this case is found to be 



77(r) ci 2F1 \^pi,-pi;qi; 
k + 2m + 2 
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2k + 2 



92 
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k — 2m 
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1 + e 



2fc + 2 



j^s 2k+2 
d) 



(16) 



These wavefunctions strongly resemble the k = case. For m > the second term is unphysical and requires one 
to set C2 — 0] likewise for to < — 1 one requires that ci = 0. The termination of the hypergeometric series requires 



Pi, 2 



Vod 
2k + 2 



= ^ + 91,2, 



(17) 



where TV is a non-negative integer which corresponds to the number of non-zero nodes in the wavefunction. There 
are no additional requirements for the regularity of the wavefunction. 



0.5 



-0.5 




(a) 



10 



r/d 



0.5 



-0.5 





1 


i 


(b) 



10 



r/d 



0.5 



-0.5 





(e) 



10 



0.5 



-0.5 







i 


(d) 



10 



r/d 



r/d 



FIG. 1: (Color online) Radial wavefunction components for the first two states (A^ = 0, 1) with angular momentum m = 1 for 
the Lorentzian dot {k = 0) and the first ring {k = 1): (a) fc = 0, iV = 0; (b) fc = 0, iV = 1; (c) k = 1, N = 0; and (d) fc = 1, 
A'^ = 1. Solid (dotted) lines correspond to components xa (ixs)- Insets: shape of the probability density for each state. 



The radial wave- function components for to = 1, the first angular momentum to permit confined modes, are plotted 
in Fig. [T] for both the Lorentzian dot and the simplest k = I ring. As an example the first fully integrable mode, which 
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corresponds to node number N 
wavefunction 



with m = 1 for the Lorentzian dot, has the following normalized two-component 
(r/d) e'" [ 1 



*(^'^) = \\\- 2 

dV [1 + (r/d)2]2 



(18) 



IV. PROBE MICROSCOPY 



The immediate physical relevance of the electrostatic Lorentzian potential becomes apparent by considering a 
standard image charge problemi^ We construct a toy model of an STM tip in which the tip is treated as a spherical 
charge held above the graphene plane at the charge neutrality point (i.e., E — Q\ depicted schematically in Fig. [2l 
which produces the potential 



C^stm('') 



eQtip 
47rK 



(19) 



where k is the permittivity of the air, Qtip is the total charge accumulated on the STM tip and hi and /i2 are the 
distances from the metallic surface to the graphene layer and STM tip, respectively. Clearly the electrical potential 
energy experienced by charge carriers in the graphene plane is independent of the radius of the spherical tip when the 
charge is fixed, as follows from Gauss' law. The charging of the tip by an applied voltage is, however, dependent on 
the specific shape of the system via the capacitive coupling between the tip and the Si layer. 
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FIG. 2: Schematic representation of a proposed experimental setup in which a charged STM tip with characteristic radius i?tip 
produces confined states in the graphene sheet either by varying the charge on the tip or by fixing the tip voltage Kip and 
moving the tip vertically with respect to the flake. The presence of bound states is exhibited by peaks in the conductivity 
whenever new states appear. 

Evidently, this form of the potential created by the tip only includes the leading-order term due to point charges. 
It is relatively simple to include higher order terms that obey the boundary condition of a fixed potential on the 
surface of the tip by positioning an infinite number of imaginary point charges inside the charged tip and its imagci^ 
However these higher order terms do not significantly affect the results. 

The Lorentzian dot potential can trivially be fitted to Eq. (jl9|) by equating the potential maxima at the origin and 
the integrals of the two potentials, yielding the matching criteria, 



Vo 



and 



hj- hj ^ fh2 + hi 
a — ; In 



nhi 



h2 - hi 



(20) 



(21) 
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FIG. 3; (Color online) Comparison of a repulsive electric potential created by a charged STM tip (blue dotted line) to the 
solvable Lorentzian potential (red solid line) as a function of scaled radial distance. The potential profiles clearly are very 
similar. Chosen values of h\ and /i2 here are arbitrary; the profiles resemble that shown here for any values so long as h2 > hi. 
Uo = VohvF is the maximum value of the potential at the origin. 

The result of the fitting procedure using Eqs. ([20]) and ([2T|) is shown in Fig. [3l It is experimentally feasible to achieve 
a probe tip with a radius of curvature i?tip ~ 20nm and to set the height of the tip above the graphene plane of 
the order of lOnm without deforming or damaging the graphene sheeti^ The thickn^ of the air gap is of the order 
of hundreds of nanometers, such that the capacitance between the tip and its image is almost independent of the 
separation distance between the tip and the graphene flake, and so it can be approximated to 



Thus we expect the values of the tip voltage at which confined states occur to be given, to a good approximation, by 



where Vtip = Qtip/Ctip. As long as the tip does not come too close to the graphene plane, the voltage needed to form 
states is only weakly affected by the distance from the graphene, and instead the radius of the tip is the dominant 
factor. This simplified model suggests a tip of radius 20 nm held at height 20 nm above the graphene plane, suspended 
over an air gap of 280nm (dimensions that are experimentally feasible)^ would create confined states at integer 
multiples of 38mV. According to Ref. [l^, the emergence of confined states will be accompanied by conductivity 
peaks in the sample maintained at zero Fermi energy by the back gate. 



We have shown, using a model class of potentials, that full confinement of the charge carriers in graphene by 
electrostatic potentials is indeed possible for zero-energy states, so that the Fermi level coincides with the Dirac 
points. This setup can circumvent the effect of Klein tunnelling providing a new way to create quantum dots. The 
exact solutions shown here, describing both model quantum dots and rings, allow one to understand more general 
properties of such systems. Namely, square-integrable states are only possible when the angular momentum quantum 
number m is non-zero, and, due to the sublattice symmetry, also when m ^ —1. Thus the fully confined states are 
zero-energy vortices, which require for their existence the critical effective strength of the confining potential (the 
product of the potential strength and characteristic size). We have analyzed the feasibility of experimentally probing 
the presence of confined states in quantum dots induced by an STM tip above the graphene surface. 

The implications of our results go beyond the discussed STM geometry and can be used for a better understanding 
of the zero-energy states in undoped graphene. Indeed a circularly symmetric potential induced by impurities should 



Ctip « 47rKi?tip- 



(22) 




(23) 



V. DISCUSSION 
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also support zero-energy vortices upon certain conditions. At first glance full confinement can only be achieved by 
fine-tuning the potential depth and width. However the natural tendency of the system towards local neutrality, 
which minimizes the total energy including the energy associated with the electrostatic field, allows one to speculate 
that the screening would adjust the potential parameters in favor of the existence of zero-energy states. Notably, this 
will create non-zero charge density at the charge neutrality point in disordered samples providing another plausible 
cause of the minimal conductivity of graphenej^ The discussed non-linear screening should also ease the observation 
of fully-confined states in STM spectroscopy, as it will adjust the strength and shape of the potential to support 
zero-energy states for an interval of tip voltages instead of a fixed voltage value. 

Our results also provide some insight into the problem of zero-energy excitons in pristine graphene which can be 
understood as electron-hole vortices with non-zero relative angular momentum of constituent particles. The condition 
of relative confinement depends on the electron-hole interaction strength, governed by the dielectric constant and the 
spread of the potential which is infiuenced by either the separation from the substrate or by screening. Therefore, our 
results might be useful for the ongoing search for the excitonic Mott transition in two-dimensional systemsi^ 

It is well-known that the zero-energy bound states are intimately related to resonant scattering of low-energy 
particles. This is reflected in the Levinson theorem,"2§ which is known to hold in the two-dimensional case .^^'^^ Our 
current work deals with zero-energy states only. Their mergence with the continuum when the potential is close to 
the critical strength and their influence on scattering and conductivity is a subject of future study, for which the 
knowledge of exact solutions for model potentials should provide a valuable tool. 
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